SU(2)xU(l) unified theory for charge, orbit and spin currents 
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Spin and charge currents in systems with Rashba or Dresselhaus spin-orbit couplings are formu- 
lated in a unified version of four-dimensional SU(2)xU(l) gauge theory, with U(l) the Maxwell field 
and SU(2) the Yang-Mills field. While the bare spin current is non-conserved, it is compensated by 
a contribution from the SU(2) gauge field, which gives rise to a spin torque in the spin transport, 
consistent with the semi-classical theory of Culcer et al. Orbit current is shown to be non-conserved 
in the presence of electromagnetic fields. Similar to the Maxwell field inducing forces on charge 
and charge current, we derive forces acting on spin and spin current induced by the Yang Mills 
fields such as the Rashba and Dresselhaus fields and the sheer strain field. The spin density and 
spin current may be considered as a source generating Yang-Mills field in certain condensed matter 
systems. 

PACS numbers: 72.25.-b, 72.10.-d, 03.65.-w 



I. INTRODUCTION 

Spintronics or spin based electronics offers opportuni- 
ties for a new generation of electronic devices for informa- 
tion process and storage Q, Q • One of the recent devel- 
opments in this field is the study of the spin Hall effect, 
which has potential applications to generate and manip- 
ulate spin polarization and spin currents by applying an 
electric field [1 H 1 H H ©, EJ E E M HE M u3- 
The intrinsic transverse spin current induced by an elec- 
tric field was predicted by Murakami et. al. 13 in the 
Luttinger Hamiltonian and by Sinova et al. ,5] in two- 
dimensional electron systems with a Rashba spin-orbit 
coupling. A resonant spin Hall conductance was also 
predicted by Shen et al. in the latter system when 
a strong perpendicular magnetic field is applied. It has 
been generally agreed by now that the spin Hall conduc- 
tivity vanishes in Rashba systems with impurities in the 
absence of a magnetic field 0, @ ■ On the experimental 
side, the coherent spin manipulation and the electrically 
induced spin accumulation have been observed |rA llOllll| . 

One of the current issues in study of spin Hall effect 
is the non-conservation of the bare spin current in sys- 
tems with spin-orbit interaction, which is unfamiliar to 
the community and requires better understanding and 
interpretation. Recently, Culcer et al. 01 have devel- 
oped a semi-classical theory of spin transport and intro- 
duced a torque dipole moment to the spin current. In 
this paper, we present a unified SU(2)xU(l) theory for 
spin and charge currents. In our theory, both the Rashba 
and Dresselhaus interactions are described by an SU(2) 
Yang-Mills field, and the total spin current contains an 
additional contribution from the field strength tensor of 
the Yang-Mills field, and is conserved. Our theory pro- 
vides a microscopic understanding of the spin torque in- 
troduced in the semi-classical spin transport theory. The 
orbital- angular- momentum current, or the orbit current 
is shown to be non-conserved in the presence of an elec- 



tromagnetic U(l) field. We also derive the forces induced 
by the Yang Mills field such as the Rashba and Dressel- 
haus fields and the sheer strain field acting on spin and 
spin current. Finally we argue that the spin density and 
spin current may be regarded as a source to generate 
Yang-Mills fields. 



II. PARTIALLY CONSERVED SPIN CURRENT 

It is well known that the Dirac equation for an electron 
moving in an external potential V (as well as the Maxwell 
vector potential A) in the non-relativistic limit up-to the 
order of 1/c 2 leads to the Hamiltonian, 



H = 



1 / e a + ^ m s x A V 
2m c 



h 2 

mX 2 



V. 



(1) 



where A 



ft 2 

4m 2 c 2 



VV, B = V x A, and e denotes the 



charge of the carriers under consideration. 

The first term in Eq. indicates the Rashba spin- 
orbit coupling {e.g., X — az), which represents the dy- 
namical momentum involving the interaction of an elec- 
tron with both the U(l) Maxwell field and SU(2) Yang- 
Mills field; the second is the Darwin term [^j, the third 
term is the Zeeman energy, and the last one is a higher 
order term [2l| . 

For an electron system confined in the x-y plane, 
the strength vector of the Rashba coupling is given by 
A = b{ E ) |22, with E the electric field along the 
z— direction and b the linear coefficient. The component 
of the SU(2) Yang-Mills field potentials in this system 
can be expressed in terms of the U(l) Maxwell fields E 
and B, hence the SU(2) Yang-Mills field strength is re- 
lated to the Maxwell- field strength and their derivatives. 
This appears to be a realistic example of Yang-Mills field 
in condensed matter system. 
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To start with general formalism, let us consider the 
Schrodinger equation for a particle moving in an external 
U(l) Maxwell field and an SU(2) Yang Mills gauge field, 

ih%-*(T, t) = fftf (r, t) 
at 

H = ^— (p — - A — V A a r a ) 2 + eA + V A%T a (2) 
2m c 

where 'J is a two-component wavefunction. Clearly, in 
comparison to Eq. Q), Aq and Aq refers to + — 
and — jH B s ■ B, respectively, but our further discussion 
are valid for the general case. In this paper, the Greek 
superscripts/subscripts stand for 0, 1, 2, 3, the Latin ones 
for 1,2,3, and repeated indices are summed over. Let 
us denote the vector potential of the Maxwell electro- 
magnetic field by = (Ap, Ai), and the vector po- 
tential of the Yang-Mills field d8] by A M = i°r°, 
with r a the generators of SU(2) Lie group. Similar 
to the field-strength tensor of the electromagnetic field 
(Maxwell field), which is defined by F M „ = d^A v — d^A^, 
the field strength tensor of the Yang-Mills field is de- 
fined by F M „ = J-^t"-, whose components are given by 
) A a — h A a — nf abc A b A c 

Let s a = (s 1 , s 2 , s 3 ) be the spin operators, then the lo- 
cal spin density is given by a a (r, t) = W(r, t)s a ^(r, t). 
Using the method similar to the derivation of the con- 
tinuity equation for the charge current, one obtains a 
continuity-like equation 

(— -r)Iox)a(T, t)+(— +r ) I i x)j i (r, t)=0 (3) 
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Yang-Mills field Ai produces a torsion r/Ai x Ji on the spin 
current, which results in an additional spin precession. 

Note that the continuity-like equation for the partially 
conserved spin current Eq. © can be written in a matrix 
form, 



(G) 



with A M = nA^t ". Here l a stand for a representation 
matrices of the generators of SO (3) Lie group, namely, 



/0 ON 

4= o o i 

V0 -1 0, 



This provides a geometric picture that the non- 
conservation of spin current implies the existence of a 
non-trivial connection characterizing parallel displace- 
ments in space-time manifold, i.e., the local intrinsic 
frame rotates from point to point. 



III. UNIFIED OPERATORS FOR CHARGE 
AND SPIN CURRENTS 

As we shall see below, it will be instructive to consider 
spin and charge in a four-dimensional intrinsic space with 
three dimensions in spin and one in charge. The conti- 
nuity equation related to the conventional charge density 
p = e^ty reads 




where the spin-current density is defined as 

J(r, t) = -SRe^^sv + vs)f (4) 

with the velocity operator v given by 

dr 1 1 p 

* = Tt=-H^ H ^^-c A -^ < 5 > 

and r\T a — s a which are just the spin operators if r? = %. 
In order to avoid ambiguities, the indices i, j, k, /i, v 
refer to the components of a vector in spatial space while 
a, b, c, a refer to those in the intrinsic space (Lie algebra 
space, or spin space). Moreover, a vector in spatial space 
is denoted by a bold face while that in intrinsic space 
is specified by an overhead arrow, e.g., a = (a 1 , a 2 , a 3 ), 
A = (A U A 2 ,A 3 ), Ai = (AlAlAf), J - (J 1 , J 2 , J 3 ) 
and J 3 = (J 3 , Jf, Jf) etc.. 

Let us discuss the physics implications of Eq. J2J. In 
the static case J = 0, Eq. reduces to dtd = r/Ao x a, 
which implies that the 0-th component of the Yang-Mills 
field Ao induces a torsion on the local spin density a. 
From Eq. (3), we can also see that the spin current is 
non-conserved even in the limit .4q = provided that 
the spatial components Ai ^ 0. This is because that the 



(7) 



where j = eKe'^vvI' is the charge current density. 
We introduce a four-dimensional velocity and a four- 
dimensional current valued on the U(l)xSU(2) gauge 
group. Let r Q = (t°, t") with r° being the identity 
matrix, the four-dimensional velocity operator is defined 
as = w°r Q , which is determined by the Heisenberg 
equation of motion, 



ih ^ 



(8) 



In the above equation, we have taken the commutation 
relation between xq = ct and the Hamiltonian to be ich. 
For the Hamiltonian J5J), we can write 



= (fio, fit) = (c, v° + <r a ), 



(9) 



with ug = 1, = 0; vf = (ft - ~Ai)/m and vf = 
—rjA^/m. Using these velocity operators, a unified four- 
dimensional "current tensor operator" including both 
charge and spin degrees of freedom, is then defined as: 



ra 



1 



(10) 
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The anti-commutator guarantees the defined tensor to be 
hermitian. To make the stipulation consistent with what 
we have discussed above, a = a with a = 1, 2, 3 refer to 
spin and a = refers to charge. Accordingly, s° = er°. 
From Ea. lJT0*j l we obtain, for \i = 0, 



fa 



• J o 



Jg = cs a 



(11) 



which recovers the usual definitions of charge and spin 
densities, 

p(r) ee - J °(r) = i^tj"* = e*%, 
c c 

cr a (r) ee i (r) = itftjo^ = tft^. (12) 
For fJ, = i, we have 



which defines the charge and spin current densities, 
ji(r) ee Jf(r) = Stetftjj°t& = eSie^tfii*, 
Jf (r) ee Jf (r) = Ke^t Jf # = afte* 1- -^ + 1°^)*. 

(13) 

The continuity-like equations © can be written in a four- 
dimensional form 

Q 

— J^M^A**, J p ] =0, (14) 

where JJ M = J^t q denotes the spin-current vector which 
is a matrix valued vector, and A M is obtained from by 
the Minkowski metric tensor g^ v — diag(l , —1, — 1, —1). 
This shows that Eq. © is an SU(2) covariant exten- 
sion of traditional continuity equation. By comparing 
the above equation with the continuity equation for the 
charge current, we see that the non- vanishing terms in 
Eq. ©, or Eq. (|14fl in systems with Rashba spin-obit 
coupling is just the spin torque in the semi-classical the- 
ory of Culcer et al. [l^] • The spin torque represents the 
spin procession caused by an external magnetic field and 
the Rashba interaction. Since the rotational symmetry 
in spin space is broken, the spin density is not conserved. 
In our theory the non-conservation of spin current in the 
presence of Yang-Mills field is due to its non- Abelian fea- 
ture. 

Since the spin current defined in Eq. |@J is non- 
conserved, it will be interesting to examine the origin 
of the non-conservation by means of constructing a con- 
served " total current" . 



IV. CONSERVATION OF A TOTAL SPIN 
CURRENT 

An electron has two important intrinsic properties: its 
charge and spin. As we have illustrated, the Rashba or 



Dresselhaus coupling describes the interaction between 
an electron and some particular SU(2) gauge fields. The 
SU(2) gauge potential was also adopted in the discus- 
sion of the quantum interference of a magnetic moment 
in magnetic fields |l9j. We are therefore motivated to 
construct a theory with the gauge field coupled to spin 
and spin current in Lagrangian formalism. 

The non-relativistic Lagrangian density of the system 
in the symmetric form is given by 



£>nr 
1 

2m 
1 



[(p - -A - ?]A a T a )«f] * ■ [p — -A — r]A a T a ] * 
1 



-F F 



•T~a -77a 



(15) 



where * = 



According to the Euler-Lagrangian 
equation, this Lagrangian density gives the same equa- 
tion of motion as Eq. J5J 



9 T 



r 1 

2m 



{p--A-r 1 A a T a y + eA + r ] A%T a ¥. (16) 



For the system with a Rashba spin-orbit coupling in some 
semiconductors, the SU(2) gauge potential is antisym- 
metric Al = —A l a , namely, 



e lab d b V. 



(17) 



However, our formalism here in terms of the Yang-Mills 
field covers more general cases, such as systems with a 
Dresselhaus coupling. 

Noether's theorem suggests that the invariance of the 
system under a continuous transformations will lead to a 
corresponding conservation quantity, hence a conserved 
total current can be defined by 



rrOL 



^eeJ" 
5A^ a M 



J 



which obeys 



g 7 a 



0, 



(18) 



(19) 



where J'" refers to the contribution from the Yang- 
Mills field and Maxwell field. We adopted a notation 
A* = (A°, A°) = (Ap, A%) in Eq. JTSJl and the explicit 
components of can be easily read out as 



J ° = e* 1 *, 

J° = eUte^ViV, 

for the charge case and 

Jf = rj8te*ti(r°iS i + fi i T a )*, 



(20) 



(21) 
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for the spin case. These results are consistent with the 
unified four-dimensional currents we have introduced in 
previous section. Eq. 120JI and Eq. 1)21(1 are just non- 
relativistic limit of ib^f^ip an d ^Iblfil 11 ^ 1 , respectively, in 

which ib —\ refers to Dirac fcrmions. 

W 

Now let us turn to the extra part J'" For a = 0, J' 
is null and the charge current j = e^te^^v^ (i.e., </?) 
is just the total current and hence is conserved, which 
is consistent with Eq. 0. Note that the field-strength 
tensor of the SU(2) Yang-Mills field J 7 " depends on both 
the derivatives of the potential A® and the potential it- 
self. This property is very different from the Maxwell 
field, which only contains the derivatives of the poten- 
tial Afj,. Thus the contribution from the Yang-Mills field 
J'u = 7/e ohc 'A h v T^ v is not null, namely, 



J'o = -r\Aj x £j, 

J'i = rjAo x Si + r)e ijk A 3 x B k , 



(22) 



where S and B are the Yang-Mills "electric" and "mag- 
netic" fields. The above density of (intrinsic) angular mo- 
mentum J' and its flow J' i cancel the non-conservation 
term brought about by the bare spin current after sub- 
stituting them into the continuity equation 1(19(1 for the 
total current. 

The Lagrangian in Eq. ((15(1 is shown suitable for an 
electron subjected to the Rashba spin-orbit interaction 
((17(1 . while that formalism can be constructed for other 
systems with spin-orbit interaction. It is our conjecture 
that the non-conservation of the spin current of the elec- 
tron system may always be associated with corresponding 
Yang-Mills field. In this picture, the Yang-Mills field in- 
duced by spin and spin current carries a fraction of the in- 
trinsic angular momentum, in analogy to the the Maxwell 
field induced by charge and charge current which car- 
ries orbital angular momentum. Finally, we argue that 
spin and spin current may be regarded as a source gen- 
erating Yang-Mills fields, which in turn interact with the 
electron. With the help of Lagrangian formalism similar 
to Eq. 1(15(1 . one should be able to study various physi- 
cal properties of the system, including the effect of the 
Yang-Mills field to the observable. Since the effect of the 
Yang-Mills field is weaker than that of the Maxwell field, 
the bare spin current has been suggested to be employed 
in detection Hal. 



V. FORCES ACTING ON SPIN AND SPIN 
CURRENT 

Using the covariant formalism, we can easily determine 
the force induced by the Yang-Mills field on the spin den- 
sity and spin current (named as spin force for simplicity) . 
Analogous to the Lorentz force evaluated by j^F^ in 
electromagnetism, the general form of the force provided 
by a Yang-Mills field is given by 

fi = J^^ = a a Si - e ijk JjBl ^ (23) 
where a a and Jf stand for the spin density and spin 
current, respectively. 

The Yang Mills fields corresponding to the Rashba and 
Dresselhaus spin-orbit couplings, the Zeeman term, and 
the sheer strain field for a system in x-y plane are given 
by: 



2^ 
V 



( B x ,B y , B z ), 



— * 2 771 

M = -jr( 0, a, 72/), 

7 2m. . 

A 2 = —{-a,-p, -jx), 



A:, 



(0, 0, 0), 



(24) 



where B = (B x ,B y ,B z ) stands for the magnetic 
(Maxwell) field, a and f3 refer to the Rashba and Dressel- 
haus coupling strengths, and 7 is related to sheer strain 
field, respectively. The "electric" Yang-Mills field is given 
by Ej = S?r a with 



^ = ^-diB + 2u n Ai x B. 
V 

The "magnetic" Yang-Mills field reads 

B 3 = -(d ia + d 2 (3y - -£(0l/3 + d 2 a)r 2 

rjh r\n 



(25) 



Am 2 , ., , -mi 
H W (P 2 -<x 2 )-^)r, (26) 



4m 



while Bi = B 2 = 0. In the above equations, a and j3 
are non-uniform in general. The force can be derived 
explicitly 



2 ^ B -a ■ d x B + [0(ff x B) t - a(a x B) 2 ] 



f 2 = ^a-d 2 B 
h = ^*-d 3 B, 



rjli 
rjh 



4*1 



J 2 K> 2 

J 2°3 



t3k?3 
J 2 °3> 



[a(a x B)i - 0(ff x B) 2 ] + j\B\ + J\B\ + J\B, 



3- 



(27) 
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The first term in these equations corresponds to the force 
due to inhomogeneity of the magnetic field, which is 
the same as in the Stern-Gerlach apparatus. Clearly, 
there are more forces acting on the spin, related to both 
the magnetic field and the Rashba and Dresselhaus cou- 
plings. Furthermore, the spin current will be subjected 
to transverse forces. Note that the inhomogeneous mag- 
netic field generates a force of Stern-Gerlach type on spin 
while the non-uniform Rashba and Dresselhaus coupling 
introduces a force on the spin current. In what follows 
we consider some special cases. 



VI. ORBIT DENSITY AND ORBIT CURRENT 

It is worthwhile to investigate the continuity-like equa- 
tion for orbit density and orbit current. We can define 
a local density of orbital angular momentum w a (r,i) — 
Ivj/t^avji _|_ I(L Q vp)t\]/. Here L is the dynamical angular 

momentum, L a = e abc Xb(p c — (e/c)A c — riA a T a ) for the 
system of an electron moving in a Maxwell field and a 
Yang-Mills field. The orbit density so defined is gauge 
covariant. Using the Shchrodinger equation J5J) we can 
derive a continuity-like equation 



Uniform Rashba and Dresselhaus fields 



d 



d_ 

dxi 



I? 



:ibc 



■ h b- r cv°i/ ~ u 



(31) 



If the Rashba and Dresselhaus coupling strengths are 
uniform (a and /3 are constants) and the magnetic field 
is absent, we simply have 



h = o. 



(28) 



Clearly, the forces arising from the Rashba and Dres- 
selhaus couplings are along the opposite direction. The 
magnitudes of the forces are related to the perpendicular 
component (a = 3) of the spin current only. Eq. (|28|l was 
also obtained in a semi-classical approach |2fi| . 



B. Non-uniform Rashba and Dresselhaus fields 

If a [25l | and (3 are non-uniform, there will be a trans- 
verse force whose magnitude is related to the in-plane 
components (a = 1, 2) of the spin current. For example, 
in the case |a| = \(3\, and B = 0, we have 



h = 

h = 
h 



1m 

— [4(8 ia + 8*0) + 4(diP + d 2 a)] , 
1m 

-—[Jlid^ + d 2 /3) + 3\{d x p + d 2 a)] , 



0. 



(29) 



Pure sheer strain field 



For electrons subjected to a sheer strain field ,;24] only, 
which corresponds toa = /3 = _B = but 7^0, the 
forces are given by 



where the flow of the orbital angular momentum, namely 
the orbit current is given by 



1 



It(r,t)=^I^--I9 



4e 



1=0 P 



(32) 



and 



I? = (v.^ + L a Vi)/2 

v { = ^ - (e/c)Ai - v A a T a )/m 



(33) 



In Eq. i(3"T|) — F ci refers to the field-strength tensor 
of the Maxwell field and to that of the Yang-Mills 
field; J° = j i refers to the charge current and J? the 
spin current. Clearly the orbit current l|32|l so defined 
is also gauge covariant. Note that the definition of the 
orbit current is not a simple extension of the spin current 
Q by replacing the spin operator by the orbital angular 
momentum operator. Actually, there is one more term 
in the orbit current, which involves change density p. 

It is worthwhile to point out that the third term in 
Eq. (|31|l is related to both charge current and spin cur- 
rent. In the absence of Maxwell and Yang-Mills fields, 
the third term in Eq. I|31l) vanishes and the orbit current 
is thus conserved. Using notations lu = (u>i, 0J2, ^3), 
Ii = {I}, if, If) and r = (xi, x 2 , £3), we can write 
Eq. I|31|) in the following form 



d d - 

— u+ — U =r x F + r x / 



(34) 



The physics implication of the right-hand side represents 
torque produced by Lorentz force F due to Maxwell field 
and the aforementioned spin force / due to Yang-Mills 
field. 



(30) 



VII. SUMMARY 



where 7 oc C^/{he) in the notation of Ref [2J|. Clearly, 
there is a transverse force acting on the spin current. 



We have introduced a four-dimensional charge and spin 
current tensor for systems coupled with Yang-Mills fields. 
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The Rashba spin-orbit interaction and Dresscnlhaus in- 
teraction can be regarded as particular Yang-Mills fields. 
The current tensor is related to the SU(2)xU(l) gauge 
potential. We have also provided a precise definition 
of orbital angular momentum current. Using the La- 
grangian formalism, wc have constructed a conserved to- 
tal current, which consists of a conventional bare spin 
current and a non-vanishing term contributed from the 
Yang-Mills field. The latter provides a microscopic in- 



terpretation of the presence of a spin precession resulted 
in the non-conservation of the bare spin current. Wc 
have derived a general formula describing the forces act- 
ing on the spin and spin current. We have proposed that 
the spin density and spin current can be regarded as a 
source generating Yang-Mills fields in a similar way as the 
Maxwell field generated by charge density and current. 

The work is supported by NSFC grant No. 10225419 
and by RGC in Hong Kong. 
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